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ji^, Motivated by the recently observed anomalous large dimuon charge asymmetry in 

pH ' neutral B decays and the unexpected large CP phase in the mixing-induced CP asymmetry 

for Bg — )• J/^p(j) decay, we study the effects of the anomalous tensor couplings to pursue 

possible solution. With the constraints from the obsevables (f)i^^'^ , a^^i and AM<j, the 

. parameter spaces are severely restricted. Numerically, we find the anomalies in B'^ — B^ 

!>■ ■ 

I mixing system could be moderated simultaneously by the contributions induced by the 

00 ' 

f-*) ■ color-singlet or color-octet tensor operators with their respective nontrivial new weak 

phase <j)Ti ~ 28.0° (15.6°) or (pxs ~ —62.1° (—74.2°) and relevant strength parameters 
gri « 6.7(8.5) x lO'^ or grs ~ 1.9(2.4) x 10^^ ^he CP-violating phase cpi^''"'' = 
X ; -0.77l[]-^?(-2.36l[j-g), respectively. 
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The flavor changing neutral current (FCNC) processes, such as 6 — )■ s transitions, arise 
only from loop effects within the Standard Model (SM), and therefore are very suitable for 
testing the SM and probing its various extensions. Among the many decay models induced by 
6 — s transition, — oscillation is of great importance, which is governed by a Schrodinger 
equation 



^ \Bsit)) ) \ 2 / ^ 

with the mass matrix M'^ and the decay matrix . The mass and the width differences between 
the light and the heavy Bs mass eigenstates {Bl and Bh) are 

AM, = Mlj-Mi = 2\M{^1 Ar, = - T?, = 2\T{^\co^ct>,, (2) 

with the phase 0, = arg(— M^2/ri2)- Another independent observable is the like-sign dimuon 
charge asymmetry for the semileptonic decay of mesons, which is defined by 

A',, ^ ^ = fi,ai + ^al, , with a,^^ = lm^ = i^sin0,, (3) 

where N^^ and A'^^ represent the numbers of events containing two mesons decaying 
semileptonically into two positive or two negative muons, respectively. 

Within the SM, the theoretical evaluations for these quantities have been fully studied [H 
[21 E] . The recent updated SM predictions for the mass and the width differences are AM, = 
(19.30 ± 6.68)ps-^ and AF, = (0.096 ± 0.039)ps-i which agree well with the CDF Col- 
laboration result for AM, [1] and Heavy Flavor Averaging Group (HFAG) averaged data for 
AF, [5] 

AM, = 17.77 ± 0.10 ±0.07ps"^ CDF (4) 
AF, = 0.154^^^7^ or -0.154l°:°5™ps-\ HFAG (5) 

respectively. While, some recent measurements of CP violating observations in Bs mixing 
system present some anomalies, which may be the hints for new physics (NP) and motivate the 
search for new source of CP violation. 

One of the hints for the existence of new CP violating source in — system is an 
unexpected large CP phase in the mixing-induced CP asymmetry for Bs — t- J / ^0 measured by 
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CDF [HI E] and DO [SI [H] Collaborations. Averaging the data from CDF and DO, two possible 



solutions ( named SI and S2 for convenience) for the CP phase 







are given by HFAG [5] 



0.77+°f, (SI), 
\ -2.36l°:g (S2) . 
However, within the SM, the CP violating phase is predicted to be small 

^W^SM) = 2arg[-y,,y,:/K,K:] ^ -0-040 , (7) 

which deviates from experimental data by more than 2.5(T. 

Another hint for new source of CP violating is the recently observed anomalously large 
CP-violating like-sign dimuon charge asymmetry in semileptonic B hadron decays reported by 
DO Collaboration [10]. Using the data corresponding to 6.1/6~^of integrated luminosity, DO 
Collaboration reports the result [10] 



= (9.57 ± 2.51(stat) ± 1.46(sys)) x 10"^ , (8) 

which differs by 3.2a from the SM prediction (— 2.3to 6) ^ 10~^ [2]. With the known values for 
mixing parameters of B^^g system [TI], DO Collaboration obtains [TU] 

All = (0-506 ± 0.043)af^ + (0.494 ± 0.043)a^^ . (9) 

Using the known experiment value a^^ = (— 4.7±4.6) x 10"'^ [5], Eq. ([9]) leads to a^i = (— 14.6± 
7.5) X 10~^ ||10j. Combining CDF Collaboration measurement of A^i [T2] and DO Collaboration 
direct measurement of a*; [12], one can get the average value 

a^, = (-12.7 ±5.0) X 10-^ (10) 

which deviates from the SM prediction (2.06 ± 0.57) x 10~^|2j by 2.5cr. If such a large a*^ is 
confirmed, it would imply the existence of a significant new source of the CP violation beyond 
the SM in the — 5^ mixing system. 

In many extension of the SM, various new couplings could be generated and represented 
by new four-quark operators in effective theory. One kind of them is the anomalous tensor 
operators, which are helpful to resolve the abnormally large transverse polarizations observed 



^Theoretically, the phase cjyi^^'^ is different from 0^ appeared in Eq. ([3|) (see Ref. [2] for detail). 



in i? — )■ (f)K* decay, as well as the large B{B r]K*) [T31 [13 [IS] and have attracted many 
attentions recently [T7j. In this Letter, motivated by the aforementioned anomalies in — 
mixing system, we shall pursue possible solutions through a set of anomalous FCNC tensor 
operators in a model independent way. 

The full set of the tensor operators responsible for the B^ — B^ mixing belonging to the 
TLL, TRR and TLR sectors read 

TLL: O^f = s,a^,(l-75)&.^,^^'^(l-75)&,, O^^ = s,a^,il - ^,)b,s,a''''{l - (11) 
TRR: O^f = s,a^,(l+75)Mia^"(l + 75)&,, O^g^ = s,a^,(l + 75)&,s,or^"(l + 75)&i ; (12) 
TLR : O^f = Sia^,{l - ^5)ks,a^''{l + 75)&i , O^f = ^^^^,(1 - ^5)bjSja^''{l + l^)hi , (13) 

where i and j are color indices. For the TLL operators, because the renormalization group (RG) 
evolution and the parametrization of the matrix elements are the same as TRR sector, its 
contributions could be easily obtained by replacing the corresponding NP parameters; As for 
the TLR operators, the hadronic matrix elements relevant to the two operators O^f g are zero 
in vacuum insertion approach. So, in this Letter, we shall pay our attention only to the TRR 
operators for simplicity. 

To begin with, at electro-weak (EW) scale, the four fermion interactions responsible for 
5° — 5° mixing induced by the tensor operators 0|?f and O^^ could be described by the 
effective Lagrangian 

^J// = -GWriO^^ + 9lsO^^) + h.c. , (14) 

where the new effective chiral h — s FCNC tensor couplings gxi and qts are generally complex 
and could be written as 

9T1 = IgTile'"^^' , 9T8 = \9Ts\e"^^' • (15) 

Using the color-singlet operators basis introduced in Ref . [18] , the effective Lagrangian Eq. (HM 
could be rewritten as 

^Iff = -Gl [{g'^, + lg'rs)0?^ - GgkO^^] + h-c. , (16) 
through the Fierz identity 

0^3« = -60ff + l0^f, (17) 
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where = Si{l + 75)6jSj(l + 75)&j. Then, including both the SM and the tensor operators 
contributions, the full expression of the effective Hamiltonian for AB = AS = 2 transition at 
rub scale could be written as 

Kffif^b) = ^.rnl^iVM' [C^'if^M'^ + C^^M^ + C|^{^^b)0|?] + h.c. , (18) 

where Oy^ = Sj7^(l — 7s)&jSj7^(l — 7s)6j, Cy^^fib) and C^^g^i^^b) are the Wilson coefficients 
at the scale ^b = ^b- In Eq- f lTS]) . the first term corresponds to the SM part, and the last two 
terms are the contributions induced by the tensor couplings. At the EW scale fiw = Mw, the 
Wilson coefficients read 

C^\fiw) = So{x,) + ^^^[S,{x,) + FMSo{x,) + B,So{xt)], m (19) 

^^^^^^^ - ^ wzr ' ^ 

^^^^^^^ - iv,bV,iY ■ 

For convenience, in the following we will absorb the factor into the effective couplings 
parameters qti and qts- The Renormalization Group (RG) evolution of these Wilson coefficients 
from the fiw scale down to fib scale have been fully developed in Ref. [TB]. With the NLO r] 
factors given by Ref. [18] , the explicit expressions for the Wilson coefficients at fib scale are [18] 



Csi^(^fe) I _ ( [^ll(/^6)]sRR [^12(^fe)]sRR I I C§^{flw) 
Criif^b) I \ [^2l(/ife)]sRR [^22(/^fe)]sRR / V C^^if^w] 



(22) 



The hadronic matrix elements corresponding to the operators in Eq. ( 1T8|) could be parameterized 

as 

(O^'-ipt)) = jml./I.Silw), (23) 



Finally, in terms of the effective Hamiltonian given in Eq. (|T8|) . the off-diagonal Mfg in the 
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Table 1: The values of the theoretical input parameters. 



as{M^) = 0.1184, Gf = 1.16637 x 10"^ GeV-^ mw = 80.399 GeV ,mB^ = 5.366 GeV, 
m,(2GeV) = 0.101+°:°^9 GeV , m,(m,) = 1.27+°;°^ GeV ,mb(mfe) = 4.20+°;^^ GeV , 
vole ^ 4 79+0.19 ^poie ^ ^ ^ 3 GeV .mArnt) = 164.8 ± 1.2 GeV , 


[II1II9] 


A = 0.804 ± 0.010 , A = 0.22535 ± 0.00065 , p = 0.135 ± 0.040 , f] = 0.374 ± 0.026 , 


m 


/b, = (0.231 ±0.015) GeV, 

= 0.86 ± 0.02+°:°^ , B2{fib) = 0.83 ± 0.02 ± 0.04 , B3{fib) = 1.03 ± 0.04 ± 0.09 . 


[SUES] 


Bs mass matrix is given by 






(26) 


where the first term is the SM contribution for and could be rewritten as [1] 




2mBM2iSM) = ^M^{VtbVZf{BBjl)ml^VBSoixt). 


(27) 


In the SM, the off-diagonal element of the decay matrix have been fully evaluated in 
Refs. [21 E] 


rusM) = -[xlm + 2A,A„r- + XlT^^] 

= -[A^r- + 2 Ai A. (r- - r-) + a^ (r- - 2r- + n^")] 


(28) 


with the CKM factors \i = V*.Vib for i = u,c,t. The explicit expressions for r^^'""''"'^ 


could 



be found in Refs. [21 E] • It is important to note that the NP operators Eq. fill I) considered in 
this Letter can significantly affect but not ri2, which is dominated by the CKM favored 
b — CCS tree- level decays within the SM. Hence = ^i2i^^) holds as a good approximation. 

With the relevant theoretical formulae given above and the values of the input parameters 
summarized in Table [H we now proceed to present our numerical analyses and discussions. 
Including all of the theoretical uncertainties, our results of the SM predictions are listed in 
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Table 2: Numerical results for AM,[ps~^], Ar,[x 10"^ pg-i], 0^/'^'^ and a^;(xl0~3) within the 
SM and the anomalous tensor couplings. 

Exp. SM Case A Case B 

SI S2 SI S2 

AM, 17.77 ±0.12 18.25 ±3.69 17.77 ±0.20 17.77 ±0.20 17.77 ± 0.20 17.77 ±0.20 

Ar^ 15.4t^;o (or -15.41^;°) 9.5 ± 4.0 5.9 ± 4.7 -4.8 ± 4.4 5.5 ± 4.3 -4.8 ± 4.5 



-5.. 

-0.77to 3? U -2.36lg jg -0.040 ±0.004 -1.08 ±0.39 -2.07 ± 0.42 -1.06 ± 0.36 -2.05 ± 0.40 
-12.7 ± 5 0.026 ± 0.009 -5.7 ±1.4 -5.7 ±1.4 -5.4 ±1.1 -5.6 ±1.3 



the third column of Table HI One may easily find that our results AM^ = 18.25 ± 3.69 ps~^ 
and Ar^ = (9.5 ± 4.0) x 10~^ps~^ agree well with the experimental data 17.77 ± 0.12ps~^ 
and (15.417;q) x 10^^ps~\ respectively. However, the predictions (pi^'^'^ = —0.040 ± 0.004 and 
a'li = (2.6±0.9) X 10^^, which agree with the former SM results (for example, the ones in Ref. [2]), 
significantly deviate from experimental results — 0.771;q;37 U —2.36~^q^I and (— 12.7± 5) x 10~^, 
respectively. In the following, we shall pursue possible solutions with the anomalous tensor 
couplings. 

For the convenience of analysis, with the central values of the theoretical input parameters 



listed in Table. [H the amplitude of — mixing could be written as 

{Bs\Kjfifit)m = A^''{SM)+Af^^ + Af^i, (29) 

Av^{SM) X 10^^ = (6.26 - iO.25) , (30) 

A^^ X 10^^ = -0.148 X \gTi X lO^pe*^'^^^ , (31) 

A^R X 10^^ = 1.83 X \gTs X lO^pe^^-^^* . (32) 



It is found that the contributions of the tensor operators Af^^ and A^^ with IsfTi.rsI ~ 

0(10-2) 

could be comparable with the SM contribution Ay^{SM) and may resolve the anomalies in the 
B'^ — Bg mixing. 

With the central values of the theoretical input parameters, the dependance of the observ- 
ables (f)i^^'^ and a^i on the new phases 0Ti,r8 with different \gTi,Ts\ values are shown in Fig. [TJ 
From Fig. [T] (a), we find that (pi^^'^ could be brought to the experimental level of SI with 
(pTi ~ 20° and \gTi\ ~ 6 x 10-^. A larger \gTi\ is demanded to enhance (p'i^^'^ to the experi- 
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Figure 1: The dependance of and a^^ on the new phases (j)Ti,T8 with different |(7ti,T8| 

values. For each plot, the irrelevant NP parameters are set to zero. The dashed lines correspond 
to the error bars (Icr). The legends for figures (b) and (d) are the same as the ones labeled in 
figures (a) and (c) respectively. 



mental result of S2 than that of SI. Meanwhile, with the same values of \gTi\ and (j)Ti, from 
Fig. [1] (b), it is interesting to note that a^i also could agree with the experimental result by the 
contribution of Al^^. For the effects of Af^^, as shown in Figs. [T](c) and (d), a similar situation 
could be found with a smaller l^f^sl ~ 2 x 10"^ and a negative phase (pTs ~ —65°. Therefore, 
the large discrepancy between the SM prediction and experimental measurement of both (pi^'^'^ 
and a^i could be moderated simultaneously by the contributions of O^^ and/or O^^ operators. 

From Figs|T](b) and (d), one may find that large \gTi\ (~ 9 x 10^^, for example) or {qtsI (~ 
3 X 10~^, for example) would lead to much larger contributions to MI2, however, they are 
disfavored by the anomaly a'li, which also can be found from Eq. ([3]). Furthermore, such large 
IqtiI and/or {qtsI are also very fragile under the constraint from the well measured AM,,, which 
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lgTll[xlO-^] I^8|[X10-^] 

m 0) 

Figure 2: The allowed regions of the tensor couplings parameters I^ti.tsI and 0ti,T8 in case A 
and B under the constraints from (pi^^'^ (green region), a^^ (pink region), AMg (blue region) 
and their combination (meshed region). 



Table 3: Fitting results for the tensor couplings parameters Is'ti.tsI [>< 10 ^] ^-nd 0Ti,r8[°]- 



NP pars. 


Case A 


Case B 




SI S2 


SI S2 


\9ti\ 


6.7 ±1.8 8.5 ±1.4 




4>T1 


28.0 ±7.1 15.6 ±7.3 








1.9 ±0.4 2.4 ±0.3 


0T8 




-62.1 ±6.9 -74.2 ±7.1 



SM prediction agrees well with the experimental data. Therefore, combining the constraints 
from (pi^^'^, all ^Mg, the NP parameter space is restricted very much. So, in order to 
evaluate the possible strength of the NP effects, we shall perform a detailed numerical analysis 
in the following. 

Due to the cancelation between the contributions of Aj^i and Aj^g, which can be found in 
Eqs. fl3T]) and fl32l) . tensor couplings qti and qts are hardly to be well bounded without any 
further simplification. In order to evaluate their respective effects, our following calculations 
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and discussions are divided into two simplified color scenarios named case A and case B, which 
correspond to with only color-singlet or -octet tensor operator considered, respectively. In each 
case, our fitting for the tensor couplings is performed with the experimental data on (pi^'^'^, 
a^; and AM^ within l.GScx (about 90% CL.) as constraints and the theoretical inputs allowed 
within their respective uncertainties listed in Table. [T] As for the observable AF^, due to its 
weak constraint on the tensor couplings with its large experimental error bar, we leave it as 
our prediction, which could be tested by the refined measurement in the forthcoming years. 
The allowed regions of the tensor couplings parameters IgriTsl and (priTs in the two cases 
under the constraints from (pi^^^, a^i, AMg and their combination are shown in Fig. [21 and 
the corresponding numerical results are listed in Table [3l where the two solutions SI and S2 
correspond to the two experimental results of (pi^'^'^ given in Eq. Our theoretical results 
for the observables are summarized in Table [2j 

Case A: only color-singlet operator 

In order to evaluate the effects of color-singlet tensor operator, we neglect the NP contributions 
involving qts- Under the constraints from (pi^^'^, a^i and AMg, Fig. [2] (a) shows the allowed 
spaces of the tensor coupling parameters. We find all of these constraint from (pi^^'^, a'li and 
AMg are improtant to get the restricted NP parameter space. Numerically, the tensor coupling 
parameters are seriously bounded to \gTi\ = (6.7 ± 1.8) x 10"^ ((8-5 ± 1.4) x 10"^) and (p^i = 
28.0° ± 7.1° (15.6° ± 7.3°) in SI (S2), which are nontrivial to moderate the discrepancies of 
(pi^'^^ and a^i between the SM predictions and the experimental measurements. 

Corresponding to the values of \gTi\ and (pTi in SI and S2, our theoretical results of the 
observables are listed in the fourth and the fifth columns of Table [2l We find our result of (pi^'^'^ 
agrees well with the measurement within la. Meanwhile, a^i is enhanced to the experimental 
level by the contributions from Af^^, but it also deviates from the experimental data by about 
l.lcr, which is easily understood through the following analysis. 

From Eqs. ([2]) and ([3]), one can obtain 

at, = ^2 sin <Ps, (33) 

for which, we can take AMg ^ 17.77 due to the well measurement and agreement with SM pre- 
diction, and |F^2l = |ri2l5M — 0.051 due to the assumption of the negligible NP contributions 
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to TI2. With the central value of the experimental result for a^y as input, using Eq. (15^ . one 

can get 

- 12.7 X 10"^ ~ 2 sin 0, sin 0, ~ -2 , (34) 

which is obviously out of the allowed sin 0^ values. So, a^i is hardly to be enhanced to around 
—12.7 X 10"'^ by NP contributions. In fact, because the CP violating phase is very small 
within the SM 0f*^ ~ 0.004 and would be dominated by the NP contributions, the relation 
^SM+NP ^ ^J/ixP ^ _o.77+o-29 y -2Mtoil would be a good approximation within NP mod- 
els P]. From Eq. fl551) . one may find that any of the various NP models without providing 
significant modification to would be impossible to give a large negative a^^ as —12.7 x 10"'^. 
While, one may note that there are large uncertainties in the experimental measurements for 
and AF^, which will hinder a definite conclusion. 

Due to the relation </)f*^+^^ ~ (f)i^^^, which leads to cos0f^^+^-^ < cos^f^^ ~ cos(0.004), 
our prediction AF^ = (5.9 ±4.7) x 10~^ps~^ in SI listed in Table [2]is smaller than the SM result 
(9.5 ± 4.0) X lO^^ps"^ and the averaged experimental data {l^A^ij'^) x 10~^ps~^ [5]. Taking 
into account of the large experimental error and theoretical uncertainties, they are still in good 
agreement. Recently, CDF and DO Collaborations have updated their measurement 

AF, = 0.075 ± 0.035 ± 0.010 ps-\ CDF [7] (35) 
AF, = 0.15 ± 0.06 ± 0.01 ps"^ DO (36) 

One may find that the CDF result is much smaller than the DO measurement, but close to 
our result (5.9 ± 4.7) x 10~^ps^^. The combination of CDF and DO measurements is not 
available until now, which will be very important for further constraining the parameter space. 
Additionally, since there are two possible intervals for ^f^'^'^: > -0.77+°;^^ (Sl)> -tt/2 and 
-it/2 > -2.36^0:29 (S2)> -tt, our prediction AF^ > in SI and AF^ < in S2 are hsted in 
Table [21 which agree with the experimental results 15.4ly;Q or —15.415 4, respectively. 

Case B: only color-octet operator 

In this case, we consider the NP contributions induced by color-octet tensor operator solely. 
Our fitting results for the tensor parameters (pTs and IqtsI are listed in the fourth and the fifth 
columns of Tableland shown in Fig. [2] (b). Due to the different sign of Af^§ and , which 

11 



could be found in Eqs. and 1^, a negative 0t8 = -62.1° ±6.9° (-74.2° ± 7.1°) in SI (S2) 
is demanded by the constraints from (pi^^'^ and a^i, which also can be seen from Figs. [1] (c) 
and (d). Meanwhile, the strength of the color-octet tensor couplings \gTs\ ~ 1.9(2.4) x 10~^ is 
much smaller than the color-singlet one \gTi\ ~ 6.7(8.5) x 10~^ in SI (S2). The situation of our 
numerical results and the effects of color-octet tensor operator for the observables in this case 
is similar to the ones of case A. 

In summary, we have studied the recently observed anomalies in — mixing with a set 
of possible tensor operators in a model-independent way. It is found that both (pi^^'^ and Ogj 
could be bridged to the experimental data by the contributions of color-singlet or color-octet 
tensor operators with new weak phase (pTi ~ 28.0° (15.6°) or (pTs ~ —62.1° (—74.2°) for the 
two possible solutions for the CP-violating phase (f)i^^'^ given by CDF and DO measurements. 
So far, it is unknown that these tensor operators with coupling and weak phase determined 
in this letter could be generated in which realistic NP models. It surely deserves further 
studies. It is noted that the error bars in the measurement of AF^, (pi^^'^ and a^^ are very 
large, which hinder us from deciphering color-singlet or color-octet tensor operator responsible 
for the anomalies. With the running LHC-b experiment, the refined measurements of these 
observables are expected to confirm or refute the NP effects. 
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